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Abstract. We investigate electron transport through a mono-atomic wire which is tunnel coupled to two 
electrodes and also to the underlying substrate. The setup is modeled by a tight-binding Hamiltonian and 
can be realized with a scanning tunnel microscope (STM). The transmission of the wire is obtained from 
the corresponding Green's function. If the wire is scanned by the contacting STM tip, the conductance 
as a function of the tip position exhibits oscillations which may change significantly upon increasing the 
number of wire atoms. Our numerical studies reveal that the conductance depends strongly on whether 
or not the substrate electrons are localized. As a further ubiquitous feature, we observe the formation of 
charge oscillations. 



PACS. 05.60.Gg Quantum transport 
Quantum wires 



73.23.-b Electronic transport in mesoscopic systems - 73.63.Nm 



1 Introduction 

Mono-atomic wires of metal atoms fabricated on a sur- 
face are the ultimately small conductors and may be used 
in nanoelectronics to connect nanodevices such as quan- 
tum gates, qubits, or nanotransistors. Thus their elec- 
tronic properties are of crucial interest. One-dimensional 
mono-atomic wires can be fabricated using mechanically 
controlled break junctions Q3H]- Such wires are freely sus- 
pended and, thus, do not interact with any substrate. For 
the same reason, they are somewhat unstable, and it con- 
sequently represents a challenge to form long wires. Simi- 
lar but more stable structures can be fabricated on vicinal 
surfaces and investigated with scanning tunneling micro- 
scopes (STM) [5 111115] ■ see Fig. [TJ Well ordered and even 
longer examples are double stranded gold wires grown on 
silicon vicinal surfaces such as Si(335) and Si(557) j3USJ|B]. 
The geometry and the electronic structure of these setups 
are sufficiently stable such that measurements can be re- 
peated many times. Notice that STM experiments mainly 
focus on the electron transport from the STM tip to the 
surface (perpendicular transport) and, thus, are not con- 
cerned with the transport from one end of the wire to the 
other end. The latter type of STM experiments would re- 
quire some modifications of the setup as is discussed in 

Refs. [zususumi • 

The conductance of ideal or disturbed wires has been 
investigated both experimentally and theoretically. For ex- 



ample, it has been predicted that the conductance of an 
atom chain depends on whether the number of atoms is 
even or odd [11II12II13II141IT5] . These even-odd oscillations 
have been confirmed experimentally pQ- Conductance os- 
cillations with larger periods may occur as well [To1[T7, 18 , 
IT51 . They stem from a Fabry- Perot like resonance of elec- 
trons with Fermi wavelength in the chain, which eventu- 




elektrode S 




Electronic address: tomasz.kwapinski@umcs.lublin.pl 



elektrode S 

Fig. 1. Schematic view of a quantum wire on a surface (elec- 
trode S) contacted at sites 1 and 4. The contacts is established 
by (a) one fixed electrode and one STM tip or (b) two STM 
tips. In both configurations, the right electrode can be moved. 
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ally changes the chain filling factor P^lll8llT5] . Moreover, 
the formation of charge waves inside a wire was observed 
for both, non-magnetic [3D] and magnetic wires [21]. The 
fact that the resonance condition depends on the presence 
of further atoms beyond the contact also leaves its finger- 
prints in the conductance, where interference effects can 
be observed [2^[23l[24l[25l [26 l [27 ] . 

In this paper we consider setups in which one lead is 
realized by a movable STM tip by which an atom of choice 
can be contacted; see Fig. [T] We describe the atom chain 
by a tight-binding Hamiltonian and obtain the conduc- 
tance within a Green's function approach [15,18,28,29, 
30,31,32,33 . It will turn out that the conductance oscil- 
lations are influenced by small leakage currents from the 
wire to the substrate. Our substrate model, the spatial 
separation of the wire atoms is taken into account by con- 
sidering local tunnel couplings. A relevant parameter of 
this model is the Fermi wavelength of the substrate elec- 
trons, which allows one to interpolate between two limits: 
In the one limit, each wire atom couples individually to 
a substrate with localized electrons, while for the other 
limit model, the substrate represents a reservoir of de- 
localized electrons. These two setups may be interpreted 
as an insulating and a conducting substrate, respectively, 
or the coupling to a molecule with according localization 
properties [51] , 

The paper is organized as follows. In Sec. [2] we present 
our model and a scattering formalism. Moreover, we derive 
an analytical expression for the retarded Green's function 
of a wire coupled to a surface. With these expressions at 
hand, we investigate in Sec. [3] both conductance oscilla- 
tions and charge oscillations. The main conclusions are 
drawn in Sec. 2J 

2 Theoretical model and formalism 

We consider the setups sketched in Fig. [1] which both con- 
sist of a quantum wire connected to two metallic elec- 
trodes. The wire may exchange electrons also with the sur- 
face which, thus, represents a further, weakly connected 
electrode. One of these electrodes may be fabricated by 
epitaxy or grown on the surface and is fixed |10j . The 
other electrode is movable and contacts an atom of choice. 
Alternatively, both electrodes may be realized by an STM 
tip as is sketched in Fig. QJb). 

The model Hamiltonian for a wire with N atom sites 
can be written in the form H = Hq + H tun , where 

N 

H o = ^2 e ka al a a ka + ^^etafai (1) 

ka=L,R,S i=l 

describes the electrons in the wire and in the leads. Elec- 
tron transitions between the leads and the wire are estab- 
lished by the tunnel Hamiltonian 

Htun — 22 ^kLO>kL a m + VhRO^R^ 
k k 

N-l (2) 

+ + n - c - 

1=1 



Here, m and n label the atom connected to the left and to 
the right lead, respectively. It is worth mentioning that if 
the left STM electrode couples to the first wire site, m = 1 
and the right electrode to the last site, n = N, the system 
corresponds to the break junction geometry of Refs. [TJ 
IT3llT^IT5"irrg] . The operators <Zj and a/" create and annihi- 
late, respectively, an electron at site i = 1, . . . , A, while 
ctka and aj£ are the according leads operators. The tunnel 
matrix elements Vki enter the expressions for the current 
only via the spectral densities F l = 27r^ fc \Vki \ 2 S{s— £ke), 
£ = L, R, which we model within a wide-band approxima- 
tion as energy independent. With the above Hamiltonian 
we assume that electron-electron interactions do not lead 
to correlation effects and can be captured by an effective 
shift of the onsite energies. Then both spin directions are 
independent of each other, such that the spin need not 
be considered explicitly. For Au or Pb chains on vicinal 
silicon surfaces, these conditions are met reasonably well. 
Generally, this should hold for non-magnetic wire atoms 
[35] . It has also been shown that electron-electron correla- 
tions do not change the period of conductance oscillations 

In order to describe electron leakage from the wire to 
the surface, we consider the surface as a further, weakly 
coupled electrode. Thus, we introduce the wire-surface 
tunneling Hamiltonian 

flw-a =J2 V ik a h a i+ h - C -> ( 3 ) 
fc 

where VA — V£ exp(ikRj) is the tunnel matrix element 
for atom j [551159] . The phase factor reflects the position of 
atom j and has the consequence that the leakage depends 
on the spatial separation of the atoms. Assuming equal 
distances a between neighboring atoms, while mainly sub- 
strate electrons with Fermi wavelength play a role, we ob- 
tain Vj k = Vfc exp(ikpja). As for the leads, the influence 
of the substrate can be subsumed in a spectral density. It 
will turn out that the relevant quantity reads 

s s sm(k F a\i-j\) 

k F a\i-j\ ' (V 

with the effective leakage strength r s . A formally similar 
coupling and spectral density has been used to describe 
decoherence of spatially separated qubits coupled to a 
bosonic environment [40,41,42 . While the fermionic case 
can still be treated within scattering theory, the bosonic 
model gives rise to memory effects which may be consid- 
ered within a non-Markovian master equation approach 

gMI5] 

Two limiting cases are worth being discussed: (i) If 
fcpd 3> 1, the spectral density is rather small unless i = j. 
This means that we can employ the approximation rfj = 

r s Sij . This describes a substrate with a very short mean 
free path such as a semi-conductor or an insulator. Then 
an electron that tunnels from a particular atom to the 
substrate can re-enter only at the same site. Obviously, 
this scheme corresponds to a model in which each wire 
atom is coupled to an individual additional electrode, (ii) 
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In the opposite limit, kpa <C 1, we obtain rfj = r s . 
Physically this means that the wire electrons tunnel to 
a delocalized substrate orbital, as is the case for metallic 
surfaces. 

In order to obtain the linear conductance between the 
electrodes L and R and the local density of states at site i, 
pi, one needs to compute the Green's function for the total 
Hamiltonian. The linear conductance at zero temperature 
is given by the Landauer formula [2^129151111511155] 

G=^T(E F ) = ^r L r R \G r mn (N,E F )\\ (5) 

where T(Ep) is the electron transmission at the Fermi en- 
ergy which we choose to be Ep = 0. Using the equation of 
motion for the retarded Green's function G r mnl one finds 
the elements G r mn from the relation G r mn = (A _1 ) m „, 
where the matrix A is given by 

(A)ij = (e — - Vi(Sij + i + S l+ ij) 

.r^ .£l (6) 

2 2 2 

Below, when presenting specific results, we will always as- 
sume that both leads couple equally strong to the wire, 
yielding E L = r R = r, and that all onsite energies and 
intra-wire tunnel matrix elements are position indepen- 
dent, £i = £o and Vi = V. These assumptions are quite 
reasonable for a wire consisting of one atom species in an 
equidistant arrangement on the surface. 

In the absence of wire-surface tunneling, i.e. for r s — 
0, the matrix Aij becomes tri-diagonal and its inverse, 
the Green's function, can be computed analytically. After 
some algebra we find 

G r mn = (-V) n - m det A»~ n { det Aq + i^ m<n 

r (7) 

_ _detA™- 1 det^- n detv4[ 1 l - m - 1 } , 

where 

£»i n2 = det Af' 1 det A^ nl + det Af' 1 det A^ n2 , 

(8) 

while Aq denotes the tri-diagonal N x N matrix for r s — 
r = 0, i.e. for the isolated wire. The determinant of this 
matrix can be evaluated to read detA^ = V N Ufq{(f)), 
where u m {4>) is the iVth Chebyshev polynomial of the sec- 
ond kind and 4> = arccos{(e — So/2V)} plays the role of a 
Bloch phase [2D]. Note that det A^ = l and det A\ = s-e . 

For a wire-surface coupling in the limit (i), the ad- 
ditional self-energy is proportional to the unit matrix, 
= r s 8ij, such that A^ 1 can still be computed ana- 
lytically. Then we obtain again the Green's function ([7]), 
but with e — eo replaced by e — £o — ir s /2. 

The local density of states at wire site i is determined 
by the retarded Green's function G r u owing to the relation 
Pi(s) = — Im (cof Aff/ir det A N ), where coi Af( denotes 
the algebraic complement of the matrix Af t , the so-called 



cofactor. The charge localized at site i can be obtained 
by integrating the local density of states up to the Fermi 
energy, 

Pi{e)de. (9) 

-oo 

Note, that due to wire-surface coupling, analytical formu- 
las for the local density of states or charge density Qi can- 
not be simplified further. Thus Eq. ([7]) is the most general 
analytical expressions for the retarded Green functions for 
arbitrary wire length. 



3 Conductance oscillations 

The conductance of a quantum wire is governed by the 
electron wave functions at the Fermi surface. In particular 
for short wires, it may even be such that a single orbital in 
the relevant energy range dominates. Its overlap with the 
electrodes may change significantly with the length of the 
wire and, thus, the conductance changes as well. Typically 
this change appears as periodic oscillation [TTllT21[T51[ni 
IT5lH6l[T7lll8lfT9] . For a break junction geometry, i.e., when 
the left electrode is connected to the first atom (m = 1), 
and the right electrode to the last atom (n = N), the 
period M of the oscillation can be determined analytically: 
Writing the transmission T(Ep) in terms of the Green's 
function (J7J), one imposes that this transmission for a wire 
of length N must be the same as for a wire of length 
N + M . This results in the condition cos(ttI/M) = (Ep — 
£o)/2V, where I = 1,...,M - 1 [H]. The same relation 
holds for a wire coupled to a surface. 



3.1 Quantum wire isolated from the substrate 

Before addressing the influence of a substrate, we study 
the conductance oscillations of a wire that couples only 
to the electrodes, but not to the substrate. In particular, 
we focus on the influence of the wire atoms beyond the 
STM tips. The presence of these atoms modifies the wave 
functions of the wire electrons and, consequently, it may 
affect the overlap between the relevant wire states and the 
electrodes. In all our numerical studies presented below, 
we use the Fermi energy Ep = as reference point and 
assume a fixed left electrode i.e. m = 1. All energies are 
measured in units of the wire-electrode coupling r, such 
that formally E L = r R = 1, while the conductance G is 
plotted in units of the conductance quantum Go = e 2 /h, 
such that it becomes identical to the dimensionless trans- 
mission T(Ep). 

Figure [2J depicts the conductance of wires with various 
lengths N as a function of the tip position n; cf. Fig.QJa). 
The onsite energies are position independent, Si — Eq, 
and chosen such that they satisfy the oscillation condition 
for period M — 2,3,4 (panels a, b, and c, respectively). 
A most prominent feature is that the amplitude of the 
emerging oscillations depends strongly on the total wire 
length N, or put differently, on the number of atoms be- 
yond the right electrode. A closer look at the results for 
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Fig. 2. The conductance (in units of [e 2 /h] so that G = T) 
as a function of the STM tip position n for the wire lengths 
N = 8, 9, 10, 11 in the absence of wire-substrate coupling, i.e., 
for F s = 0. The onsite energies are (a) eo = 0, (b) eo = V, 
and (c) eo = v2V. These values correspond to conductance 
oscillations with periods M = 2,3, 4. The intra-wire tunnel 
matrix elements is V = 4F. The lines serve as a guide to the 
eye. 



periodicity M = 2 (panel a) reveals that the amplitude 
changes with each additional wire atom from a large to 
a small value and back. For odd N, the conductance os- 
cillates with a very large amplitude, whereas for an even 
number of atoms, these oscillations are very small. Note, 
that the conductance obtained for any even N is hardly 
distinguishable form the one for N = 10. The same holds 
for any odd N and ./V — 9. 

For the periods M — 3 and M — 4, we find large oscil- 
lation amplitudes for the wire lengths N = 5, 8, 11, . . . and 
N = 3,7, 11, . . ., respectively. For other lengths, the con- 
ductance still oscillates, but its value never exceeds 10% 
of the conductance quantum. The bottom line of these 
numerical investigations is that we observe strong oscil- 
lations with period M provided that the wire has length 
N = Mk — 1, where = 1,2,3, ... is any natural number 
[T5] . Thus for periods M = 2, 3, 4 as considered in Fig. [3J 
maximal amplitudes of the conductance oscillations are 




N=9 



N=10 



i=2 i=3 



Fig. 3. Local density of states at atom wire sites i = 1, 2, 3, 4 
for wires that consist of (a) N = 9 and (b) iV = 10 atoms. The 
onsite energies are eo = 0, such that the conductance obeys 
oscillations with period M — 2, i.e. even-odd oscillations. The 
other parameters are as in Fig. [5] The horizontal lines mark 
the Fermi energy. 



observed for 



N = 




for M = 2, 
for M = 3, 
for M = 4. 



(10) 



Some insight about the physical origin of the con- 
ductance oscillations is provided by considering the local 
density of states Pi{s). For a case with even-odd oscilla- 
tions (M — 2), this is shown in Fig. [3] for the wire sites 
i = 1,2,3,4. The magnitude of the conductance can be 
understood from the fact that for N = 9 (panel a), the 
density of states possesses a peak at the Fermi energy. 
Thus, electrons from the Fermi surface of the left elec- 
trode can tunnel to the wire resonantly, which facilitates 
transport. For length N = 10 (panel b), by contrast, the 
density of states at the first site vanishes, such that trans- 
port is practically blocked. The density of states at the 
further sites changes with period M = 2. Applying the 
same arguments to the right STM tip explains that the 
conductance must oscillate with the same period. In the 
second case (even N), however, the value of the conduc- 
tance is small due to the small local density of states at 
the Fermi energy on the left STM tip . This directly trans- 
lates to a small oscillation amplitude of the conductance. 
In the same way one can explain the conductance oscil- 
lations with other periods. In general, the conductance is 
maximal when both STM electrodes are connected to sites 
that possess a large local density of states at the Fermi 
level. 

A possible application of our results is the experimen- 
tal estimate of the onsite energy £o which strongly influ- 
ences the oscillation periods. However, since also the wire 
length N and the tip position n have significant impact on 
the conductance, it would be necessary to perform several 
measurement with wires that differ only in length but are 
identical otherwise. 
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Fig. 4. The conductance in units [e 2 /h] as a function of the 
STM position n for a wire consisting of TV = 9 atoms with 
onsite energies eo = 0. The wire-surface coupling strength is 
r s = (dashed line), Is = 0.2r (thin solid), and F (solid). 
The other parameters are as in Fig. [2] (a) Short wavelength 
limit kpa — oo, such that each atom couples to an individual 

(b) fea = 0, such that all atoms 



environment and Zy = F s 8 l3 



couple collectively to the substrate and, thus, rh = _T" 



3.2 Influence of the substrate 

The influence of the wire-substrate tunneling can be ap- 
preciated in Fig. U where we compare the conductance 
oscillations for two values of r s with those obtained in 
the absence of the substrate, i.e. for r s — 0. Let us first 
consider the limiting cases in which the substrate electrons 
are perfectly delocalized (kpa = 0) or perfectly localized 
(kpa = oo). We find two significant features: The oscil- 
lation period is practically not influenced by the leakage, 
while the oscillation amplitude decreases with increasing 
coupling strength. 

Let us now turn to the more realistic intermediate case 
of finite /cjro and its influence on the conductance oscil- 
lations, Fig. [5] Typically this parameter is in the range 
k,F<i = 1 ... 10 [39]. The onsite energies are chosen such 
that the oscillation period is Al = 2 or Al = 3. The re- 
sults interpolate those for the limiting cases discussed in 
the previous paragraph. Thus, we can conclude that the 
value of leaves the oscillation decay qualitatively un- 
changed, despite the fact that the quantitative difference 
may be significant. 

In Fig. [BJ we compare the conductance decay for per- 
fectly localized substrate electrons (kpa = oo) and per- 
fectly delocalized electrons (kpa — 0) with the one for the 
intermediate value kpa = 1.3. One notices that the de- 
cay is weaker the more localized the substrate electrons 
are. It is even such that for kpa = 0, the conductance 
does not decay entirely, but converges to a finite value. 
This reminds one to the incomplete decay of entanglement 
and coherence between delocalized qubits coupled to sub- 
strate phonons 41,42 . A qualitative explanation for the 
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n(N=9) n(N=11) 

Fig. 5. The conductance in units [e 2 /h] as a function of the 
STM position n, for a wire with N — 9 atoms and eo = (left 
panels) and for TV = 11 atoms and Eo = V (right panels) for 
intermediate wavelength such that kpa — 5 (upper panels) and 
kpa — 1 (lower panels). All other parameters are as in Fig. 0] 
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Fig. 6. Decay of the conductance oscillation with the wire- 
substrate coupling strength F s for a substrate with delocalized 
electrons (kpa — oo), with localized electrons (kpa = 0) and 
for kpa — 1.3. The other parameters are TV = 9, n = 5, and 
£o = 0. 



observed dependence on kpa is that an electron that is lost 
to a delocalized substrate orbital may tunnel back to its 
former state at any site. In the case of a substrate with lo- 
calized electrons, by contrast, the lost electron may tunnel 
back only to the very same wire site. The rates at which 
these processes occur should differ roughly by a factor TV, 
i.e., by the number of wire atoms. This crude estimate 
agrees roughly with our numerical results provided that 
r L and r R are of the same order. 



3.3 Charge oscillations 

We already mentioned that the conductance oscillations 
relate to oscillations of the charge density, which can be 
observed with STM techniques [46 . Similar charge waves 
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Fig. 7. Localized onsite charge Qi as a function of the atom 
position i for a wire consisting of N — 18 atoms with onsite 
energies (a) £q — V and (b) eo = V%V. The substrate Fermi 
wavelength is kpa = oo (thin solid line), kpa = (thick solid 
line), and = 5 (dashed). The dotted lines are obtained in 
the absence of the substrate (_T S = 0), while the other curves 
are the substrate-wire coupling is r s — F. The other parame- 
ters are as in Fig. [21 

have been predicted for a wire with break junction ge- 
ometry 20 2Tj . Thus, one may suspect that the leakage 
affects the charge oscillations as well. Figure [7] depicts the 
localized onsite charge as a function of the wire atom po- 
sition i for different surfaces and the oscillation periods 
Al = 3 (panel a) and M — 4 (panel b). Interestingly, 
with increasing leakage, the charge oscillations fade away 
stronger than the conductance oscillations. Also here, the 
influence of the leakage is more pronounced the more local- 
ized the substrate electrons are. A further common feature 
is that tunneling to the substrate affects only the oscilla- 
tion amplitude, while the period remains the same. 



4 Conclusions 

Using a scattering approach, we have studied conductance 
and charge oscillations of a quantum wire in contact with a 
movable electrode. As a particular feature, we considered 
electron leakage from the wire to various types of sub- 
strates. Our model for the latter allows for both strongly 
localized, weakly localized, and perfectly localized sub- 
strate electrons. 

As a main feature, we have found that both the con- 
ductance oscillations and the localized onsite charge are 
fading away due to the influence of the substrate. Inter- 
estingly, this influence is weaker the more localized the 
substrate electrons are. For our model, the localization 
parameter kpa, i.e., Fermi wavelength times the distance 



of the wire atoms leads to a monotonic transition between 
the limiting cases of perfect localization and delocaliza- 
tion. In all cases, the oscillation amplitude is smaller than 
in the absence of the substrate. Thus, leakage represents 
an obstacle for the experimental observation of conduc- 
tance oscillations. Nevertheless, a piece of good news is 
that leakage does not influence the oscillation period. This 
is rather encouraging, because it implies that conductance 
oscillations can be observed with wires that are grown on 
surfaces, despite their unavoidable contact to the wire. 
We thus are confident that our results will stimulate STM 
experiments with wires grown on vicinal surfaces. 
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PH acknowledge support by the DFG via SPP 1243, the excel- 
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German-Israeli Foundation (GIF, grant no. I 865-43.5/2005 ). 
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